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Abstract Solar active region (AR) 11283 is a very magnetically complex region and it 
has produced many eruptions. However, there exists a non-eruptive filament in the plage 
region just next to an eruptive one in the AR, which gives us an opportunity to perform 
a comparison analysis of these two filaments. The coronal magnetic field extrapolated 
using a CESE-MHD-NLFFF code ( |Jiang & Feng|2013] l reveals that two magnetic flux 
ropes (MFRs) exist in the same extrapolation box supporting these two filaments, re¬ 
spectively. Analysis of the magnetic field shows that the eruptive MFR contains a bald- 
patch separatrix surface (BPSS) co-spatial very well with a pre-eruptive EUV sigmoid, 
which is consistent with the BPSS model for coronal sigmoids. The magnetic dips of 
the non-eruptive MFRs match Hct observation of the non-eruptive filament strikingly 
well, which strongly supports the MFR-dip model for filaments. Compared with the non- 
eruptive MFR/filament (with a length of about 200 Mm), the eruptive MFR/filament is 
much smaller (with a length of about 20 Mm), but it contains most of the magnetic free 
energy in the extrapolation box and holds a much higher free energy density than the non- 
eruptive one. Both the MFRs are weakly twisted and cannot trigger kink instability. The 
AR eruptive MFR is unstable because its axis reaches above a critical height for torus in¬ 
stability, at which the overlying closed arcades can no longer confine the MFR stably. On 
the contrary, the quiescent MFR is very firmly held by its overlying field, as its axis apex 
is far below the torus-instability threshold height. Overall, this comparison investigation 
supports that MFR can exist prior to eruption and the ideal MHD instability can trigger 
MFR eruption. 
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1 INTRODUCTION 

As a leading cause of space weather, coronal mass ejections (CMEs) are closely correlated with solar 
filament eruptions. A recent statistic study shows that more than 70 percent of filaments eventually erupt 
and result in CMEs ( [McCauley et al.|2015| >. Therefore, it is of great importance for space weather fore¬ 
casting to understand why the filaments erupt and predict when they are likely to erupt. To answer these 
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questions, one needs the information of a key parameter, the invisible coronal magnetic field behind the 
filaments, which is believed to play the primary role in supporting the filaments and characterizing their 
stability. 

Accordingly, a variety of theoretical models have been proposed to explain the initiation of filament 
eruptions with the coronal magnetic disruption as the basis (see, e.g.,[Forbes et al .|2006) [Aulanier et al. 
|2010[|Schmieder et a 1. 1 201 3j|Aulanier|2014| and references therein). The physical mechanisms causing 
the eruption can be classified as, e.g., the breakout model (Antiochos et al. 1999| >, the tether cutting 
model (Moore et al.|200l]i as well as ideal MHD instabilities (e.g., kink instability and torus instability, 


|Hood & Priest|1981||Velli et al.|1990||Torok et al.|2004l|Tordk & Kliem|2005HKliem & Tdrok|20061 > 

and the non-ideal MHD instabilities (e.g. tearing mode instability given by Wu et al.|2000|i. 

It is understood that the pre-eruption, stressed core field in the corona is usually kept to be stable 
by an envelope/overlying field anchored strongly at the photosphere, thus, the essences of all eruption 
initiation models are to find a way how the balance between the outward magnetic pressure of the core 


field and the downward magnetic tension of the overlying field can be broken. Hence, Antiochos et al. 


( 1999) > proposed the breakout model basing on a quadrupolar magnetic field configuration, and that the 
eruption is triggered by reconnection at a coronal null above the sheared core arcade, which removes the 


overlying flux and then allows the core to escape, similar to a streamer and flux rope model (Wu & Guo 
1997; Wu et al.|1997 1 . The tether cutting model usually refers to a complete process of the formation 


and eruption of a magnetic flux rope (MFR) from an initial simple, sheared bipolar magnetic arcade 
(Moore & Roumeliotis 1992; Moore et al. 20011. In a strongly sheared bipolar field, current sheet forms 
above the photospheric polarity inversion line (PIL) and results in tether-cutting reconnection between 
the arcade field lines, which progressively transforms the arcade core into a rising twisted MFR. The 
reconnection is further forced by converging motions toward the PIL, flux diffusion and cancellation 
in the photosphere ( |Amari et al.||2003a|b[ Linker et al. 2003 ). Such processes continuously add flux 
to the MFR and in turn weaken the photospheric anchorage below the MFR, until at some instant, the 
remaining arcade tethers are too weak to prevent the eruption of the MFR. In this case, it is different 
from the breakout model since the magnetic reconnection works at low altitude below the stressed field. 
The ideal MHD instabilities mainly include the kink instability for a MFR which occurs, if the twist, 
a measure of the number of windings of the field lines around the rope axis, exceeds a critical value, 
leading to a helical deformation of the MFR’s axis (Torok & Kliem 2005; Torok etal.|2010| l and the torus 
instability (|Khem & Torok 2006), which is equivalent to a catastrophic loss of equilibrium (Forbes 1990; 
Forbes & Isenberg 1991 Kliem et al.|[2014), occurring when the MFR reaches an unstable threshold 


determined by decay speed of the envelope/overlying magnetic field that stabilizes the MFR. On the 
other hand, the non-ideal MHD instability was performed by Wu et al. ( |2000] l, who have investigated 
the formation of observed plasma blobs (She eley et al.|1997) > with multi-polarity magnetic field topology 
due to tearing mode instability. 

Of all the mentioned theoretical models, it appears that only the ideal MHD-instability one can 
provide a set of specific parameters, i.e., the twist degree of the MFR and/or the decay index of the 
envelope field, that quantifies the condition when the magnetic configuration might erupt. Thus this is 
a very attracting advantage that will promisingly make itself most popular in the community of space 
weather research. We note that such model is built up on a basic element, the MFR, which is assumed to 
exist prior to the eruption. This is indeed supported by several lines of evidence. The coronal sigmoids 
often appearing in CME-productive ARs ( |Rust & Kumar|1996 Hudson et al. 1998 Canfield et al.|1999) l 
indicate sheared and twisted magnetic configurations, thus MFR, prior to the eruption; Furthermore, the 
MFR can explain well the supporting of filaments because the dip structure in MFR provides a reservoir 
where the cool dense material can be held against gravity (Mackay et al. 2010||Guo et al. 2010b; Su et al. 


20111, and thus long-lasting quiescent filament strongly suggests a correspondingly lasting MFR; On the 


other hand, with different numerical techniques and force-free models, many authors have successfully 
reconstructed from vector magnetogram MFRs that are consistent with the observation signatures like 
sigmoids and filaments (Canou & Amari|2010 Cheng et al.|2010 Guo et alT||20 10b Jing et ak|2010 


|Guo et al.|2013[ |Jiang et al.|2014a Jiang et al. 2014b| Ruan et al.|2014 Amari et al. 2014[ ); moreover, 

direct observations of pre-existence of MFR have also been reported frequently in the SDO era with the 
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help of high-quality multi-wavelength data of the AIA (Cheng et al.|2012[|Zhang et al.|2012| |Li & Zhang 


2013a|b|c| |Patsourakos et al.|2013| |Yang et al.|2014| |Chen et al.|2014). Furthermore, the ideal MHD- 

instability theory has been extensively applied in the investigations of filament eruptions (jRomano et al. 


2003 

Williams et al. 2005 Green et al. 20071 Schrijver et al. 2008 

JLiu 2008 

Guo et al. 2010a 

Cheng 

et al. 

2011; Nindos et al. 2012] 

Xu et al. 2012f|Liu et al. 2012, Cl 

leng et al. 

[2013 Kliem et a 

. 2013 


events are in good consistent with the theory, especially the theoretical threshold values for instability 
are matched strikingly well in the data-constrained/data-driven modeling results (e.g.,|Kliem et al.|2013f 
|Guo et al.|20fOal | Jiang et al.|2014a||2013) . 

In this paper we provide further evidence for the pre-existence of MFR and the ideal MHD- 
instability mechanism for the eruption of MFR by a comparison study of an eruptive filament in 
AR 11283 and a nearby non-eruptive filament. Both filaments have been studied independently in our 
previous work (Jiang et al,|2014a||Jiang et al.|2014b) l using coronal field reconstructions, while here by 
comparing them with each other, we provide further insight in the different nature between the eruptive 
and non-eruptive filaments. We find that both of the filaments are supported by coronal MFR. The MFR 
corresponding to the AR eruptive filament is much smaller compared with the non-eruptive MFR, but 
contains significantly more magnetic free energy than the non-eruptive one. Both MFRs are weakly 
twisted and cannot trigger kink instability. The eruptive MFR is unstable because its axis reaches above 
a critical height for torus instability, while the non-eruptive MFR is very firmly held by its overlying 
field, as its axis apex is far below the torus-instability threshold height. 



Fig-1 Left: a large AIA-304 image at 22:00 UT on September 6, showing the location of 
AR 11283. Middle: closer view of the region. The top-right small panels show the eruptive 
process of the AR filaments with three time snapshots, and bottom-right small panels show 
the BBSO-Ho observations of the non-eruptive, large filament, with the arrows point to the 
filament main body as denoted in the AIA-304 image. 


2 OBSERVATIONS AND CORONAL MAGNETIC FIELD MODEL 


AR 11283 produced two X-class flares (an X2.1 on September 6, 2011 and an X1.8 on September 7), 


which have motivated many studies (Wang et al. 


2012 Petrie 2012 

Zharkov et al. 

2013 

Feng et al. 

Liu et al.j2014 

shen et al.||2014 

Xu et al. 

2014 


Zhang et al.|2015 1. Here we investigate the coronal configuration on September 6 before the X2.1 flare. 


The SDO/AIA observed that an S-shaped small-scale filament (with length of about 20 Mm) formed in 
the core of AR and erupted at the onset of the flare at 22:12 UT, with impulsive rising of the filament 
matter, as shown in Figure [T] On the other hand, in the northwest of the AR, there is a large filament, 
the main body of which seen in Ha shows a slightly inverse-S shape with a length up to 200 Mm. The 
large filament exists stably for many days, and survives during the eruption of the small one, although 
























































































4 


C.-W. Jiang et al. 


its main body is very close to the flare site. The right-bearing filament barbs as seen in the Ha indicate 
that the non-emptive filament should be related with a left-handed twisted MFR, which is contrary to 
the eruptive one. 


To investigate the coronal magnetic field that supports the filaments, we have modeled the pre¬ 
flare magnetic field using our CESE-MHD-NLFFF code with input of SDO/HMI vector magnetogram 
(see also Jiang et al.||2014b). The code is based on the MHD relaxation approach and implemented 
by an advanced conservation-element/solution-element (CESE) space-time scheme ( ] Jiang et al.|2010| l. 
We solve a set of zero-beta simplified MHD equations with a fictitious frictional force to control the 
relaxation process. The code is in line with our 3D MHD Data-driven Active Region Evolution model 
( |Wu et alT||2006[ |Jiang et aT7|| 201 2[ Jiang et ak|2013| |Wu et al7||2013[ >, which has been well tested for 
reproducing both the slow quasi-static evolution and extremely fast emption process in the large-scale 
corona. The present model is performed in a Cartesian box with the z axis along the normal vector of 
the solar surface. 


3 RESULTS AND ANALYSIS 


3.1 Comparison of coronal magnetic field model with observations 


The coronal field model presents two MFRs of opposite helicity in the same extrapolation box: a small 
one of right-hand twist (see a close-up view in Figure [2]» corresponding to the eruptive AR filament, 
and a large one of left-hand twist corresponding to the non-eruptive intermediate filament. To validate 
the model result, we carried out a detailed comparison of the magnetic field and its derived features 
with the observations of the filament related stmctures. In particular, for the large-scale filament, the 
Ha observation is compared with the magnetic dips in its MFR, since there is a good Ha image of 
the filament spine and barbs and these structures are believed to consist of cold matter collected in the 
magnetic dips. While for the small AR filament, since Ha has no adequate resolution of it, comparison 
of its MFR field lines with a EUV sigmoid that is closely co-spatial with the AR filament provides a 
better way for our validation of the modeling. 

In the left panels of Figure [2j we compare the magnetic dips with the large Ha filament. The mag¬ 
netic dips are visualized by showing part of the iso-surface ( B z = 0) with B • Vf?~ >0 (i.e., locations 
where the field lines are concave up). The dips are also pseudo-colored by the height value for a better 
inspection of the long extended dip structure, since there is a large amount of localized fragmentation of 
dips (shown in deep blue) that is very close to and on the photosphere (i.e., bald patches). Evidently, the 
long extended dip reaches above 30 Mm, exhibiting an inverse S-shaped, reproducing the main body of 
the filament. Especially, there are small-scale magnetic dips emanating from the spine, and they match 
the filament barbs strikingly well. The results are in line with and provide an evidence for the MFR-dip 
model for filaments and their barbs. 


The right panels of Figure[2]show the comparison of the magnetic field lines of the AR MFR with the 
sigmoid observed in AIA-94 (temperature of 6.3 MK). The sigmoid has a thin and enhanced forward-S 
shape, indicating a right-hand twist. We find that this MFR has a bald patch separatrix surface (BPSS), 
as shown by the thick white lines in the figure, which is the MFR’s outmost surface that touches the 
photosphere. The BP is a part of the PIL on the photosphere where the transverse field crosses from the 
negative to the positive polarity (note that the BP is directly derived from the vector magnetogram), and 
the BPSS consists of all the field lines that pass through the BP. Theoretical study suggests that strong 
current sheet can form in the BPSS of MFR and produces enhanced heating along the surface, which 
manifests itself as the X-ray or EUV sigmoid (Titov & Demoulin |1999| . The most important result is 
that the observed sigmoid coincides with the BPSS field lines, but not with those near the rope axis (the 
thick colored lines in the figure). Such fact not only supports the BPSS current sheet model for sigmoid 
brightening, but also confirms the validation of our model. 
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Fig. 2 Left: Comparison of the magnetic dips of the large MFR with the large filament seen in 
BBSO Her (adapted from Jiang et al. ( |2014b) >). The top is the magnetic field lines, the middle 
is the magnetic dips (pseudo-colored by the height from 0 to 60 Mm), and the bottom is the 
Her image. The dashed box denotes region of the main body of the filament. The arrow ’S’ 
denotes a separation of the filament. Right: Comparison of the BPSS of the AR MFR with the 
AIA-94 sigmoid (adapted from Jiang et al. ( 2014a| >). The background color images show the 
photospheric magnetic flux distribution (white for 500 G and blue for —500 G), with the PIL 
denoted by the curved line. The BPSS consists of the field lines plotted as white thick rods, 
which graze the photospheric surface at the BP. The colored field lines are sampled near the 
MFR axis. The bottom are the AIA-94 image of the sigmoid with the same field-of-view as 
the top ones. In the last panel the AIA image is overlaid by the BPSS field lines. 
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Table 1 Comparison of magnetic energy contents for the two MFRs. 



AR Eruptive MFR 

Non-eruptive MFR 

Total energy 

3.67 x 10“ erg 

1.52 x 10 32 erg 

Potential energy 

2.63 x 10 32 erg 

1.28 x 10 32 erg 

Free energy 

1.04 x 10 32 erg 

0.24 x 10 32 erg 

Ratio of Free energy to potential energy 

40% 

19% 

Field volume 

4.66 x 10 4 Mm 3 

2.98 x 10 6 Mm 3 

Mean free energy density 

2230 erg cm -3 

8 erg cm -3 



Fig. 3 Free energy density integrated along the z axis. The contour lines represent B z of 
—500 G (pink) and 500 G (cyan). The small box in the large field of view shown in the 
right panel is denoted for the AR core, i.e., for the separated volume of the eruptive MFR, 
which is enlarged in the left panel. Since the energy is much stronger in the AR core than in 
intermediate region, different scales of the energy distribution are used in the two panels. 


3.2 Magnetic energy distributions 


In Figure[3]and Table 1 we compare the magnetic energy distributions of the two filament-related MFRs. 
Figure |3]shows the free energy density integrated along the 0 axis, namely 


Efree(x,y ) = dxdy 


d2 d2 

& -^pot 

87T 


dz. 


(1) 


As can be seen, compared with the intermediate region, the free energy is significantly stronger in the 
sunspot regions, i.e., the AR core, because of the strong magnetic flux there. For the large MFR, there is 
also clearly a channel of enhanced free energy along the corresponding filament, with the central section 
possessing the strongest. Table 1 gives a comparison of different energy contents and mean free-energy 
density of the two MFRs. While the two MFRs are involved in a whole magnetic configuration, they can 


be regarded as independent of each other according to the magnetic topology. As shown by Jiang et al. 


( |2014a| >, the AR eruptive MFR is formed with the emerging of a new positive sunspot into the large 
negative flux region, and a magnetic null-point fan separatrix surface (like a dome) exists between the 
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pre-existing system and the new emerged one. The AR MFR with its overlying flux system is localized 
below the fan surface with a height of about 30 Mm. As such, we roughly separate the two MFRs in the 
full extrapolation volume using a small box with height of 30 Mm for the AR eruptive one as denoted in 
Figure[3] and the rest volume for the non-eruptive one. Then the different energy contents are computed 
using these separated volumes for the two MFRs. Notably, although in spatial size the eruptive MFR is 
much smaller than the non-eruptive one, it contains most of the magnetic free energy. Consequently the 
free energy density of the eruptive MFR is far higher than that of the non-eruptive one. The ratios of 
free energy to potential energy shows that the non-potentiality of the AR MFR field is also significantly 
higher than that of the non-eruptive one. These comparisons indicates a higher possibility of disruption 
of the AR MFR than the other one. 



Fig. 4 Analysis of the stability of the two MFRs: left panels for the non-eruptive one (adapted 
from Jiang et al. (|2014b|>) and right panels for the eruptive one (adapted from Jiang et al. 


(2014a|). The top panels show the twist degrees of the MFRs, which are calculated for several 
sampled field lines around the rope axis (which is represented by the white line for the large 
MFR in the left panel and yellow line for the small MFR in the right panel, respectively). 
The bottom panels show the results for decay index. The central cross sections for the two 
MFRs are shown respectively, and these cross sections are vertically sliced along the black 
lines shown in the top panels. The streamlines in these cross sections show the 2D field-line 
tracing on the slices, which forms helical lines centered at the axis of the MFR. The arrowed 
thick lines started from the bottom and directed through the helix center (i.e., the rope axis) 
represent the paths along which the decay index is computed. Since the axis of the large MFR 
is very long, the decay index is calculated at four more points along its axis, which are denoted 
by the colored circles in the left panels. 


3.3 Stabilities of the Flux Ropes 

To further explain why the AR filament erupts immediately whereas the large filament can keep stable, 
we study the stabilities of the MFRs with the 3D coronal field in the context of ideal MHD instabilities. 
Theoretical models and simulations suggest that a coronal MFR (i.e., MFR with footpoints line-tied on 
the photosphere) confined by an overlying potential arcade is subject to two kinds of ideal instabilities. 
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i.e., kink instability and toms instability. Kink instability occurs if the twist degree of the MFR, denoted 
by T n , which measures the number of windings of the field lines around the rope axis, exceeds a critical 
value. This threshold value is estimated to be roughly 1.5 ~ 2, according to several studies (Fan & 


Gibson|2003 Torok et al.|2004||Torok & Kliem|2005| l. The torus instability occurs when the outward 


expansion of a MFR due to its “hoop force” can no longer be confined by the overlying field (also 
referred to as the external field) if the external field decreases sufficiently fast. A decay index is defined 
as 

<9(ln B e ) 


i(R) = - 


0(ln R) 


( 2 ) 


to characterize how fast the external field B e decreases with distance R from the photosphere, and it 
has been found that if the apex of the rope axis reaches a location with the decay index greater than 
a threshold n c , the MFR system is torus unstable. This value n c is expected to lie in the range of 


1.1 < n c < 2 from a series of investigations ( 


1978 

Kliem & Torok 

2006 

Torok & Kliem 

2007 Fan & Gibson 2007; Schrijver et al. 2008f 

Demoulin & Aulanier 2010 

Au 

anier et al. 12010). One 


should bear in mind that these threshold values for instabilities may change depending on many factors, 
e.g., the specific shape of the flux rope and the magnetic environments. Thus they should be used with 
cautions in the realistic cases. 


The twist degree T n of a force-free MFR can be simply quantified by integral the force-free factor 
a = J • B /B 2 along a given field line, namely. 


T n = 


47T 


adl 


(3) 


(Berger & Prior 2006 Inoue et al.|201 1 2013[ >. In the upper panels of Figure|4] we show the results for 
several field lines around the rope axis. Since in the realistic case the MFR is actually a bundle of field 
lines winding around each other, it is difficult to locate precisely the axis of the rope like in an analytic 
or idealized model. We thus approximately find the rope axis by first making a central cross section 
of the MFR, which is a vertical slice cutting roughly through the middle of the MFR in a direction 
perpendicular to the photospheric central PIL, and then assuming the rope axis as the field line passing 
through the center of the helical shapes formed by the poloidal flux of the rope (as shown in bottom 
panels of Figure [4]). The results show that the twist degrees of both MFRs are very close to one turn, 
although different field lines of them have slightly different values of twist. Such value of twist is below 
the thresholds of kink instability that have been reported. Besides we have not observed a clear rotation 
or writhe of the erupting AR filament as it rises, which would otherwise occur in the eruption of kinked 
MFR. We concluded that neither of the filaments can trigger kink instability. 

Several points are needed to be clarified with respect to study the torus instability based on com¬ 
puting the decay index. First, since it is difficult to separate the external field from the total field in the 
model, we use the potential field that corresponds to the same magnetogram as an approximation of the 
external field, following Fan & Gibson ( |2007i ) and Aulanier et al. ( [2010] ). Secondly, as mentioned, the 
axis apex of the MFR is hard to precisely determined, and we thus approximately locate it at the center 
of the helical shapes formed by the poloidal flux of the rope in the central cross section. For the large 
filament of which the axis is rather long, we compute the decay index at four more points along its axis. 
Third, since the MFR and overlying field configuration is inclined with respect to the radial (or vertical) 
direction, it will be more reasonable to compute the decay index along a path approximately following 
the inclination of the system than along the vertical direction. Last, because the field parallel to the 
rising direction actually does not contribute to the inward confining force, the decay index is computed 
for only the perpendicular component of the potential field (Cheng et al. 2 01 1[ Nindos et al.||2012) . 
Accordingly, we compute the decay index for the two MFRs and compare them in the bottom panels of 
Figure [4] As can be seen, the non-eruptive MFR locates at heights with decay index smaller than 0.5, 
which is far below the threshold for the torus instability, thus it is very firmly held by its overlying field. 
On the other hand, the AR MFR reaches slightly above the height with decay index of 1.5, which is the 


threshold of torus instability found by Torok & Kliem (2007 1 and Aulanier et al. (2010) using idealized 
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numerical simulations. Thus our result supports the threshold value as 1.5, and suggests that the quick 
eruption of the AR MFR is due to the torus instability. 


4 CONCLUSIONS 


We have performed a comparative analysis of an eruptive filament in AR 11283 and a nearby non- 
eruptive filament in the intermediate region using a zero-beta MHD-relaxation model for near force- 
free extrapolation (Jiang & Fen g|2013] l. The extrapolated coronal magnetic field based on the measured 
vector magnetogram given by SDO/HMI shows that the two MFRs support the two filaments. Validation 
of the modeling magnetic field is performed by comparing it with EUV and Ha observations, which 
shows that the eruptive MFR contains a BPSS coinciding spatially with the pre-eruptive sigmoid, and the 
non-eruptive MFR has magnetic dips matching the shape of the non-eruptive filament. By comparison 
the MFRs with each other, the following physical characteristics are found: (i) The eruptive filament is 
only one tenth of the non-eruptive filament in size, (ii) The amount of free energy of the eruptive filament 
is five times that of the non-eruptive one (iii) The mean free energy density of the eruptive filament is 
more than two orders larger than the non-eruptive one (iv) Both the MFRs are weakly twisted, thus could 
not trigger the kink instability (v) Evaluating the decay index shows that the axis of the eruptive MFR 
(filament) reaches above a critical height for torus instability, where the non-eruptive MFR is firmly held 
by is overlying field, as its axis apex is far below the threshold height. We suggest the measured mean 
free-energy density may be a good characteristic indicator for a filament eruption in addition to other 
properties. In summary, the energy storage and the trigger mechanism are both important to filament 
eruptions, and it is supported that MFR can exist prior to eruption and the torus instability can trigger 
MFR eruption. 
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